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in dense neutral three-flavor quark matter 

Kenji Fukushima^ 

'RIKEN BNL Research Center, Brookhaven National Laboratory, Upton, New York 11973, USA 

We calculate the Debye and Meissner masses and investigate chromomagnetic instability associ- 
ated with the gapless color superconducting phase changing the strange quark mass Ms and the 
temperature T. Based on the analytical study, we develop a computational procedure to derive the 
screening masses numerically from curvatures of the thermodynamic potential. When the temper- 
ature is zero, from our numerical results for the Meissner masses, we find that instability occurs 
for Ai and A2 gluons entirely in the gapless color-flavor locked (gCFL) phase, while the Meissner 
masses are real for A4,, A5, Ae, and A7 until Ms exceeds a certain value that is larger than the gCFL 
onset. We then handle mixing between color-diagonal gluons A3, As, and photon A.y, and clarify 
that, among three eigenvalues of the mass squared matrix, one remains positive, one is always zero 
because of an unbroken U(1)q symmetry, and one exhibits chromomagnetic instability in the gCFL 
region. We also examine the temperature effects that bring modifications into the Meissner masses. 
The instability found at large Ms for A4,, A5, Aa, and A7 persists at finite T into the u-quark color 
superconducting (uSC) phase which has u-d and s-u but no d-s quark pairing and also into the 
two-flavor color superconducting (2SC) phase characterized by u-d quark pairing only. The Ai and 
A2 instability also goes into the uSC phase, but the 2SC phase has no instability for Ai, A2, and 
A3. We map the unstable region for each gluon onto the phase diagram as a function of Ms and T. 

PACS numbers: 12.38.-t, 12. 38. Aw 



I. INTRODUCTION 

The rich phase structure of matter at high baryon den- 
sity has attracted much interest over decades theoreti- 
cally and phenomenologically. At sufficiently high den- 
sity and low temperature, wherever quarks feel an attrac- 
tive force as is suggested from the one gluon exchange be- 
tween quarks that are antisymmetric in color, the Fermi 
surface of quark matter is unstable against forming a 
Cooper pair whose condensation leads to color supercon- 
ductivity T|. It is well established that the color-flavor 
locked (CFL) phase is a ground state of three-flavor quark 
matter in the asymptotically high density region where 
the strange quark mass is negligible as compared with 
the quark chemical potential fi 

If it is realized in a bulk system like the cores of com- 
pact stellar objects, quark matter must be neutral in elec- 
tric and color charges ^ong as Mg/fi is negli- 
gibly small, three-flavor quark matter, which is composed 
of the equal number of m, d, and s quarks, satisfies electric 
and color neutrality on its own regardless of whether it is 
in the normal or CFL phase. When Mg comes to have a 
substantial effect suppressing s quarks, there arise vari- 
ous possibilities in forming Cooper pairs and in achieving 
neutrality, that brings intricate subtleties into the phase 
structure especially in the intermediate density region 
where Mg can compete fj,. 

When the system is normal quark matter, a finite elec- 
tron density is required to neutralize the system elec- 
trically because the number of s quarks is reduced by 
Ms ^ 0. In CFL quark matter, on the other hand, neu- 
trality can be fulfilled even without electrons rigidly at 
zero temperature Q and approximately at low temper- 
atures below an insulator-metal crossover H S. This 



is because the BCS ansatz for Cooper pairing between 
cross-species of quarks enforces equality in the numbers 
of red, green, blue, and u, d, s quarks and thus neutral- 
ity. At zero temperature, another candidate, that is, the 
two-flavor color superconducting (2SC) phase, is known 
to cost a larger free energy than the CFL phase under the 
neutrality conditions, or, it could exist for strong cou- 
pling between quarks in some density windows 0, 13,13 ■ 

Electronless and neutral CFL quark matter lasts as 
long as the Fermi energy mismatch when two quarks were 
not paired, which is given by as explained later, 

is less than the energy gap A. There should appear a 
new state of matter once the energy gain by releasing 
the pressure between cross-species of quarks surpasses 
the gap energy, which occurs when > A. Then, 

some of quark energy dispersion relations pass across zero 
and their pairing is disrupted in the corresponding mo- 
mentum (blocking) region that is from one zero to an- 
other zero of the quark excitation energy. The new state 
is called the gapless CFL (gCFL) phase 0] and, if any 
other phase transitions such as the chiral phase transi- 
tion and the transition to a crystalline color superconduc- 
tor To, 11] lie at lower densities than the gapless onset, 
the gCFL phase must show up next as the density goes 
down from the CFL side. 

The gapless or breached pairing superconductivity was 
first discussed in a non-relativistic model and was rec- 
ognized as an unstable state and its analogue in 
three-flavor quark matter was considered in Ref. (13i] . It 
was discovered later on that stable gapless superconduc- 
tivity could be possible under the constraint of particle 
number fixing in a non-relativistic model 14] (see also 
Ref. for more discussions on stability) and of elec- 
tric and color neutrality in the gapless regime of the 2SC 
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(g2SC) phase [Tl. 

Besides gapless superconductivity, there is another 
competing pos sibihty, i.e. the mixed phase, to reahze neu- 
trahty [l7l llSl fl^. and if the mixed phase has a lower 
free energy, the gapless phase would not come to real- 
ization. According to Ref. [23| it can be the case in fact 
for the g2SC phase which should be taken over energeti- 
cally by the mixed phase where normal and color super- 
conducting phases coexist, and consequently the g2SC 
phase may be less of a reality. It must deserve further 
investigation, however, especially with the screening ef- 
fects in the mixed phase taken into account to declare 
the existence or non-existence of the (g)2SC phase. In 
the case of the gCFL phase, on the other hand, the free 
energy comparison indicates that the gCFL phase should 
be stable energetically apart from a minor exception of 
the CFL-2SC mixed phase possibility, which will be pre- 
sumably ruled out, however, once the surface tension and 
the Coulomb energy corrections are included Hence, 
so far, the gCFL phase remains as a likely candidate for 
the ground state at moderate density and is considered 
to be relevant to neutron star physics 0|. 

However, chromomagnetic instability first found in the 
g2SC phase "2^ and later confirmed also in the gCFL 
phase 24^, 25} implies that the gapless phase is not sta- 
ble against perturbation of transverse (chromomagnetic) 
gluon fields and the true ground state must need some- 
thing further. It is negative Meissner masses squared, 
i.e., i mag inary Meissner masses, that the authors of 
Refs. |23l 12^ ]25j have actually revealed in the gapless 
phases. Since the Meissner mass is the screening mass 
for transverse gluons, the instability signifies spontaneous 
generation of the expectation value of gauge fields. Keep- 
ing in mind a simple relation in superconductors between 
the gauge fields and currents in the London gauge, we can 
regard chromomagnetic instability as color current gen- 
eration. (An interpretation of instability as spontaneous 
baryon current generation has been argued in Ref. |2fi|.) 
To put it another way, it is possible to interpret it as 
instability towards a state characterized by diquark con- 
densates with color-dependent oscillation in coordinate 
space 0, because gauge fields in the quark propagator 
generally reside in the form of the covariant derivative. 
This state of quark matter is a sort of crystalline color su- 
perconductors under the single colored plain wave ansatz. 
It should be worth emphasizing that the difference is only 
in intuitive pictures, and in effect, these above interpre- 
tations are equivalent, which is mathematically described 
by the gauge transformation. Still, anyway, it is contro- 
versial what the true ground state should be that super- 
sedes the homogeneous gCFL phase. 

The purpose of this paper is to compute the Debye 
and Meissner masses in the (g)CFL phase and to present 
detailed analyses on the nature of chromomagnetic in- 
stability which occurs in connection with the presence 
of gapless quarks. Our central results are condensed in 
the unstable regions on the phase diagram presented in 
Figs. HI El and El In both the CFL and gCFL phases 
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FIG. 1: Enclosed with solid thick lines is the unstable region 
for j4i,2 on the phase diagram (see Fig. 1 in Ref. and expla- 
nations therein) . The critical end-point is marked by a square 
below which the phase transition from the uSC to the 2SC 
phase or unpaired quark matter (UQM) is first order. A tri- 
angle indicates the gCFL onset located at Mg /fi — 47.1 MeV. 
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FIG. 2: Enclosed with solid and dashed thick lines are the 
unstable regions for ^4,5 and Ag,? respectively. 



there is no mixing with other gluons for Ai, A2, A4, 
A5, Ag, and ^7, among which two gluons of {Ai^A2)^ 
(A4, ^5), and [Aq.A'j) respectively have the same screen- 
ing masses. We shall write Ai^2 for instance to mean 
whichever of Ai or A2 equivalently from now on. As for 
A3, Ag, and the electromagnetic field A^, on the other 
hand, mixing between them arises and thus we have to 
refer to eigenvalues of the 3x3 mass squared matrix to 
investigate instability. It should be noted that the color 
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which turns out to be related to the g2SC instabihty. At 
smah temperatures and large M'^/fi the Meissner masses 
squared for A^^q and Aq_y are negative until the system 
reaches the phase boundary from which the system enters 
the phase of unpaired quark matter (UQM). 
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FIG. 3: Enclosed with solid thick lines is the unstable region 
for one or two of eigenmodes consisting of A3, As, and A-,. 



indices are labeled according to the conventional GcU- 
Mann matrices in color space. 

Figurenshows the unstable region enclosed with thick 
lines inside which the Meissner mass squared for Ai^2 
is negative. The phase diagram underlaid is what has 
been revealed in Ref. Q for the case of weak coupling 
yielding Aq = 25 MeV at = T = 0. The gCFL 
phase starts appearing at M"^ / ^ = 47.1 MeV which is 
well close to the kinematical estimate (Af^//i)c ~ 2A. In 
the vicinity of the phase boundary line on which one of 
three gaps vanishes, i.e., A2 — 0, there comes out a rather 
complicated structure which is not robust but strongly 
depends on a choice of the coupling strength. We will not 
present results for other coupling cases as done in Ref. [1| , 
for our purpose here is not to disclose the phase structure 
but to argue characteristic features of chromomagnetic 
instability and to settle where it becomes relevant on the 
phase diagram. 

We observe that instability grows for Ai^2 as soon as 
the gCFL state occurs at zero temperature. The insta- 
bility penetrates from the CFL region where all u-d, d-s, 
and s-u quarks make Cooper pairs into the uSC region 
where u-d and s-u quarks remain to pair, but never en- 
ters the 2SC region where only the u-d quark pairing gap 
is nonvanishing. Actually in the 2SC phase, A12 become 
irrelevant to the Higgs- Anderson mechanism, and there- 
fore the Meissner mass for Ai^2 must be zero there. 

The unstable regions for A4.5 and Aq^y are depicted 
in Fig. 13 In the 2SC phase there is no discrimination 
between ^4^5 and Aq^y and the difference is manifested 
only when quark matter lies in the CFL or uSC phase. 
Our mapping is consistent with all of the known results 
both in the 2SC phase ^ and in the CFL phase [2i |: 
no instability takes place near the gCFL onset for ^4^5 
nor Aqj, while the large M"^ / ji regions exhibit instability 



It seems to be somewhat tricky to understand the un- 
stable regions for mixed modes of A3, A^, and A^ shown 
in Fig. 13 At zero or extremely low temperatures whose 
energy scale is determined by the electron contribution as 
discussed later, the whole gCFL phase is unstable and our 
results qualitatively agree with Ref. j2^. At finite tem- 
perature there are two distinct regions where instability 
remains. The instability near the gCFL onset eventually 
disappears as the temperature grows, while the instabil- 
ity at larger M'^/fi goes further into the 2SC phase. In 
fact in the 2SC phase, A3 decouples from others and its 
Meissner mass is reduced to zero, and one of the two 
eigenmodes composed of Ag and Aj exhibits chromo- 
magnetic instability as is consistent with the findings in 
the g2SC phase 23] . For even stronger coupling when 
the 2SC (not necessarily g2SC) phase is possible at zero 
temperature (see Fig. 17 in Ref. we have numerically 
checked that instability at larger / (x starts exactly at 
the g2SC onset. 

The exact correspondence between the gapless and in- 
stability onsets is apparent only at zero temperature since 
the finite-temperature effects allow for thermal quark 
excitations which are not clearly distinguishable from 
quarks in the blocking region. An interesting manifes- 
tation of this is the unstable region in the uSC phase in 
which there are no gapless quarks below T ^ ^ MeV (see 
the dotted curve in the uSC region shown in Fig. 1 in 
Ref. 0). In a sense, at finite temperature, we can say 
that the instability penetrates into gapped sides, as is 
also observed from Fig. 5 in Ref. 25] . 

We are explaining how we have come by these insta- 
bility mappings in a numerical way in later discussions, 
starting with the following subsections in which we shall 
make a brief review of the gCFL phase and of the deriva- 
tion of the Debye and Meissner screening masses of gauge 
fields. The readers who are already familiar with these 
basics can skip most of them and jump to Sec. ^ 



A. Gapless CFL phase 



Our strategy to approach the Debye and Meissner 
masses is based on the thermodynamic potential that has 
been formulated within an effective model. In the follow- 
ing subsubsections, we discuss the model and approxima- 
tions to describe the gCFL phase, and then illustrate the 
quark excitation energies as a function of the momentum 
(i.e. the dispersion relations). 
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1. Model and approximations 

In this paper we adopt essentially the same model and 
approximations as used in Refs. [3, . The only differ- 
ence is that the Mg effect is incorporated as an effective 
chemical potential shift. This approximation is necessary 
to relate the Meissner mass to the potential curvature 
with respect to gluon source fields in a simple way. The 
model employed here is the Nambu-Jona-Lasinio (NJL) 
model with four-fermion interaction. In our model study 
we assume that the predominant diquark condensate is 
antisymmetric in Dirac indices, antisymmetric in color, 
and thus antisymmetric in flavor; 

(CC^75^-) - Aie'^"e,,i + A2e''^2^,,^ + A3e'^^3^,,3 , (1) 

where and (a, 5) represent the flavor indices {u,d,s) 
and the color triplet indices (red, green, blue) respectively. 
The gap parameters Ai, A2, and A3 describe a 9 x 9 
matrix in color-flavor space that takes the form. 
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(2) 

in the basis (rw, gd, bs, gu, rd, rs, hu, bd, gs). In the same 
way as in Refs. p, SIj we ignore diquark condensates 
which are symmetric in color, though we know that they 
do not break any new symmetry and consequently can 
take a finite value. This approximation is motivated by 
the fact that the QCD interaction is repulsive between 
quarks in this channel and such diquark condensates have 
been actually known to be insignificant quantitatively • 
In the mean-field approximation where the conden- 
sates are relevant, we only have to consider the four- 
fermion interaction in the diquark-diquark channel that 
can be generally in hand after appropriate Fierz trans- 
formation. That is, 

£ = {^{i^ + fi-f" - M)V + Ci^t , (3) 

where 

Ant = ^m^^e'^'-^e^jvCi^J") 

Here the mass matrix M is unity in color and 
diag(0, 0, Ms) in flavor {u,d,s) space in our approxima- 
tion and fj, the matrix of quark chemical potentials in 
color-flavor space. Bold symbols generally denote matri- 
ces in color-flavor space. The charge conjugation matrix 
is C = 17^7°. 



The effect of nonzero M on the particle dispersion 
relation becomes apparent in the vicinity of the Fermi 
surface, where it can be well approximated by a shift 
in quark chemical potentials by MM^ /2fj,. We shall 
make use of this prescription that actually turns out 
to be essential to allow us to relate the Meissner mass 
to the potential curvature. In the present work, as in 
Refs. _8, 9], M is not current but constituent quark mass 
and is treated as an input parameter. We know that this 
approximation provides a correct phase structure unless 
the chiral phase transition cuts deeply into the gCFL 
region under the choice of strong coupling in the chiral 
sector (28| . 

Chemical potentials are fixed by equilibration and neu- 
trality. Stable quark matter should be a color singlet 
as a whole and neutral in the electric charge under /3- 
equilibrium. Color singletness is a more stringent condi- 
tion than neutrality. It has been shown, however, that 
no energy cost is needed to project a neutral state onto 
a singlet state in color in the thermodynamic limit [2^ . 
Hence, it is sufficient to impose global neutrality with 
respect to the electric and color charges. To describe 
that, we shall consider /ie, /^s, and /is that are coupled 
to negative Q = diag(|, — ^, — -j) in flavor {u,d,s) space 
so that a positive fie corresponds to the electron den- 
sity, T3 = diag(i,-^,0), and ^Tg = diag(i,i,-|) in 
color (red, green, blue) space, respectively. Then, 9 diag- 
onal components of the effective chemical potential ma- 
trix (including the shift by MM^ /2^) in color-flavor 
space are explicitly 

211 

^ru = M — 3 Me + 2M3 + 3^87 
l^gd = + -jMe — 5/-t3 + 

_ I 1 2 AfJ 
Mbs — ^+3^6— 3M8 2^, 

_ 2 111 

Aigti — — 3^6 — 2^^^ ' 3/^87 

_ I 1 I 1 I 1 Ml 
^rs — + 3^6 + 2/^3 + 3A^8 2^1 

2 2 

l^bu — A* ~ 3^6 — 3M81 

1 2 
fJ'bd = M + 3^6 — 3M8 ■ 

_ I 1 1,1 
fJ-gs — 3 Me - 2^3 + 3M8 " 

The gap parameters corresponding to Q hi this model 
are precisely defined as 

A, - iG(V',f"Cz75e"''''e.,,V'f )• (6) 

With all these definitions, if we choose Nambu-Gor'kov 
basis as vE'(p) = {ip{p),C-ip'^ {—p))'^ , we can express the 
mean-field (inverse) quark propagator in a simple form, 

^S~\p)^(^ + ^'f J^'^^o] (7) 
\ Z75 A ^- fi-f^ J 
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after rearranging the Dirac matrices. The quark propaga- 
tor is a 72 X 72 matrix in color, flavor, spin, and Nambu- 
Gor'kov space. From zeros of the inverse propagator, we 
can read 72 energy dispersion relations Si{p). The ther- 
modynamic potential f2 is thus written in terms of e^'s 
as 



" = -8^/'*!''i;{|£i|+2Tln(l + e-l"IO} 



J = l 



(8) 



Me 



127r2 6 180 

with the electron contribution in the last three terms. 
In order to fix three gap parameters and three chemical 
potentials at each M^, fi, and T, we will simultaneously 
solve three gap equations, 

and three neutrality conditions, 

dfi dn on 



d^e dfl3 dfls 



= 0. 



(10) 



Here A in (jHJ is the ultraviolet cut-off parameter and 
we use A = 800 MeV as in Refs. la, 9]. The coupHng 
constant G is chosen to yield Aq = 25 MeV when both 
Ms and T are zero. In the present work we take the 
quark chemical potential /i — 500 MeV that roughly cor- 
responds to 10 times the normal nuclear density in the 
model. 

When Ms is not so large as to disrupt any Cooper 
pair, /Ze, ^i3, and /xg satisfying electric and color neutral- 
ity have been analyzed at zero temperature in a model- 
independent way 4ij. Two of three can be fixed as a 
function of the third that we will choose /ig here, then 



Ai3 = ^J■e 



Ms 



.Ml 
2^ 



Me 

2 



(11) 
(12) 



The important feature of the CFL phase we should note 
is that the thermodynamic potential © at T = with 
the electron contribution disregarded is independent of 



Mq 



I (Me 



M3 



that is the chemical potential for Q — Q — T3 



(13) 



V3 



It is known that Q is a generator of the "rotated electro- 
magnetism" that is never broken by any condensate of 
the form . At zero temperature Q-charged excitations 
in CFL quark matter are all gapped and the CFL phase 
is a Q-insulator 6]. 

In the presence of the electron contribution to the ther- 
modynamic potential, in reality, there is a gentle curva- 
ture on the plateau of the thermodynamic potential pro- 
vided by the last three terms in (jS]) which select /ig = 0, 
meaning zero electron density. Although it makes only 
slight changes in energy, the existence of the electron ef- 
fects is crucial in understanding the Ai^2 instability. 



2. Gapless dispersion relations 

Since the gap matrix is block-diagonal, the 72 x 
72 quark propagator matrix is block-diagonal as well in 
color-flavor space and can be divided into four parts; one 
24 X 24 part for ru-gd-hs quark pairing with Ai, A2, and 
A3 and three 16 x 16 parts for quark pairing of bd-gs 
with Ai, rs-bu with A2, and gu-rd with A3. It might 
be possible but quite hard to handle the 24 x 24 part 
analytically. Fortunately, as we will see shortly, this 24 x 
24 intricate part has little to do with chromomagnetic 
instability, so we will limit our discussion here to the 
propagator for bd-gs, rs-bu, and gu-rd pairings. 

Leaving the explicit expression of the quark propaga- 
tor until Sec. m let us elucidate here the quark energy 
dispersion relations obtained from zeros of the inverse 
propagator. In general for the 16 x 16 part involving two 
species A and B quarks the energy dispersion relation 
takes the form. 



A-B = \ [p - ^ia-bY + A 



AB 



A-B, 



£A-B = yiP + I^'-A-bY + ^AB + ^MA-S 



(14) 
(15) 



for quasi-particle and antiparticle excitations respec- 
tively, where JIa-b = 5(ma + Ms) and 5^,a-b = \{^^A - 
fis) and Aab is the energy gap for A-B pairing. In our 
definition S^jla-b = —5^b-A- Gapless dispersion relations 
come about once 



\SfiA-B \ > Aab- 



(16) 



Then, using effective chemical potentials (0) with the 
known CFL solutions Hll|l and (|12|l . we have 



^f-^gu-rd 2 Me 2 M3 Mi 

^l^rs-bu 2Me ' 4 
S^bd-gs = -jMS — ^MS 



(17) 



1 I 1 ^'t: /-ION 

M3 + 2M8 - ^ - Me - ^ , (18) 



2' 

Ml 
it. 



4m 

Ml. 

2pL ' 



(19) 



and 



f^gu-rd — f^bd-gs — f-^rs-bu 



= M-^ (20) 



6^ 



in the CFL phase. If it were not for the electron terms in 
the thermodynamic potential, fXe could lie anywhere as 
far as it does not disrupt the pairing of gu-rd nor rs-bu 
quarks which are Q-charged. In the region between the 
dashed and dot-dashed curves shown in Fig. 01 the system 
remains to be a Q-insulator, meaning the bandgap. The 
CFL solution ceases to exist at the point where these two 
curves meet, before which a first-order phase transition 
to unpaired quark matter takes place as indicated by the 
dotted vertical line (see discussions in Ref. for details) . 

By substituting = 0, as is favored in the CFL phase 
with electrons, for the above expressions, we readily re- 
alize that the bd-gs and rs-bu dispersion relations are 
identical supposing Ai = A2 — A. In fact, the bd-gs and 
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FIG. 4: Boundaries for the gapless conditions HIGH for bd- 
gs (solid), rs-bu (dashed), and gu-rd (dot-dashed) quarks. 
Because bd-gs quarks are neutral in Q-charge, the system can 
go through the solid line, though the bandgap is edged with 
the dashed and dot-dashed lines and quark matter is a Q- 
insulator between these two boundaries. With the inclusion 
of the electron contributions to the thermodynamic potential, 
/ie is zero until M^/fj, reaches the bd-gs line, and then fie 
goes up slightly below the rs-bu dashed boundary. This is 
because the Q-charged electrons must be canceled by a small 
disruption of rs-bu quark pairing. 

rs-bu quark pairings are breached at the same time when 

\SHbd-gs\ = \Sflrs-bu\ > A , (21) 

which corresponds to the fact that in Fig. 0] the sohd 
and dashed Hnes cross at /ig = 0. This coincidence of 
the gapless onset locations for bd-gs and rs-bu quarks is 
responsible for instability for Ai^2 in the gCFL phase, as 
only Ai^2 can excite a pair of rs and gs quarks at once 
(see Fig. E)). In other words, if the electron terms are 
absent and fie is chosen to be positive, say ^ 10 MeV, 
the Ai^2 instability would disappear, while there would 
remain the instability for other gluons. 

Once quark matter enters the gCFL phase with larger 
M^//i than the onset, the bd-gs and rs-bu quark disper- 
sion relations are both gapless but no longer degener- 
ated. The blocking momentum region for bd-gs quarks 
becomes wider with increasing M'^/fi as shown by the 
solid curve in Fig. \S[ and Ai decreases accordingly be- 
cause quarks within the blocking region cannot partici- 
pate in pairing. Since rs and bu quarks have nonzero Q- 
charge, the blocking momentum region for rs-bu quarks 
is not allowed to be arbitrarily wider and is determined 
by the requirement to cancel Q-charge brought by elec- 
trons whose density is ^ fj,^. That means, the width of 
the blocking momentum region for rs-bu quarks is esti- 
mated as 5p ~ fi^/Sp,^ <C /i 0. The actual form of the 



FIG. 5: Energy dispersion relations for bd-gs (solid) and rs-bu 
(dashed) quarks at M^/fi = 80 MeV. The dotted curve is one 
that has the lowest energy among the 24 ru-gd-bs dispersion 
relations. 

rs-bu dispersion relation is, as a result, kept to be almost 
quadratic with a tiny blocking region anywhere in the 
gCFL phase as seen by the dashed curve in Fig. |31 

Although the ru-gd-bs part is hard to anticipate a pri- 
ori, there is no gapless mode involved in this sector. The 
dotted curve in Fig.^lrepresents one of energy dispersion 
relations for ru-gd-bs quarks that has the lowest energy 
among them. It is clear from the figure that ru-gd-bs 
quarks are all gapped at / fi = 80 MeV and it is also 
the case for larger M^/fi. Hence, the three (solid, dashed, 
and dot-dashed) lines drawn for bd-gs, rs-bu, and gu-rd 
quarks respectively in Fig. ^suffice for judging if gapless 
quarks are present or not. 



B. Debye and Meissner screening masses 

In a medium at finite temperature and density the 
gauge fields are screened by the polarization of charged 
thermal excitations. The screening effect on the gauge 
fields differs depending on the longitudinal and trans- 
verse directions due to the lack of Lorentz invariance. 
The screening mass, which is the inverse of the screen- 
ing length, in the longitudinal direction is called the 
(chromo)clcctric or Debye screening mass, and in (color) 
superconductors in particular the transverse screening 
mass, the (chromo)magnetic or Meissner screening mass 
in QED (QCD). In the normal phase any perturbative 
calculation leads to vanishing magnetic mass in the static 
limit. In the superconducting phase, on the other hand, 
a finite magnetic screening mass results from the Higgs- 
Anderson mechanism and embodies the Meissner effect 
so that superconductors can exclude the magnetic field. 

The Debye and Meissner masses are calculated from 
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the eigenvalues of the mass squared matrix in color de- 
fined by the self-energies; 



limn^"^(^ 



= i lira (6ij - %%)U'lJu; = 0,q), 



(22) 
(23) 



where (a,/3) represent the color octet (or photon) indices 
and we defined q — q/\q], and 11^^ is the self-energy with 
the Lorentz indices (/x,j^) for the gluons (a, /? < 8) or pho- 
ton (7 instead of a, j3) or their mixing. (We always use 
the subscript 7 for the sake of meaning photon.) At the 
one-loop level at high density, the quark loop contribu- 
tion to n^^, which is proportional to fj? , is predominant, 
so we shall consider only the polarization tensor coming 
from one-loop diagrams of quarks. 

In the symmetric CFL state with Ms = 0, the mass 
squared matrix is diagonal in color and all eight gluons 
have the same screening mass. The Debye and Meissner 
masses in a color superconductor have been calculated di- 
agrammatically in Ref. f30| and we shall quote the results 
here; 



2 21 -8 In 2 2 
™D = -Q 



1 



mj^j = -m 



D 



0.859mi 



(24) 
(25) 



for all eight gluons. Here = g fi /Gn and g is the 
strong coupling constant. Note that our definition of 
is different by a factor Nf — 3 from the definition given 
in Ref. js^. 

Once mixing between gluons and photon, which ex- 
ists among Ag and in the CFL phase at Mg — 0, 
and generally among A3, As, and A^ for Ms ^ 0, is 
taken into account, the screening masses are to be modi- 
fied. They actually have been evaluated analytically only 
when Ms = 31] and after diagonalization two eigen- 
values composed of Ag and Aj read 



m 



D,88 



(^+3^) 



4e2 
3? 
4e2 



™2 



-D,77 



™2 



''M,77 



0, 



(26) 

(27) 
(28) 



level, the gluon or photon self-energies consist of the 
particle-particle, particle-antiparticle, and antiparticle- 
antiparticle excitations. (In the case of the normal phase 
the particle-particle excitations in the vector channels 
would be to be articulated as the particle-hole excita- 
tions.) The last contribution from only the antiparticle 
excitations is just negligible so that we shall always omit 
it in our discussion. 

In the 2x2 space of Nambu-Gor'kov doubling, the 
quark propagator has not only the diagonal (normal) 
components but also the off-diagonal (abnormal) compo- 
nents proportional to the gap A that connect two quasi- 
particles. 

The Debye mass arises only from the particle-particle 
contribution. It can be further divided into the Nambu- 
Gor'kov diagonal and off-diagonal parts as follows; 



["^cldiag 



where 



[ml,]otf = 



9 2 



- 1 



-ln2|m-. 



(29) 

(30) 
(31) 



The sum of these two parts properly reproduces 124|) as 
it should. 

In evaluating the Meissner mass, on the other hand, 
it is essentially important to note that the particle- 
antiparticle excitation produces a significant contribution 
comparable to the particle-particle one. The Meissner 
mass consists of the diagonal and off-diagonal parts. 



that is further split into the particle-particle (p-p) 
particle-antiparticle (p-a) parts. 



["iM]diag(p-p) 



2['Tl^j]diag(p-a) 



with 



['™M]diag(p-p) — ol^^cldiag , 



["^M]diag(p-a) = T;'^1 : 



(32) 

and 
(33) 

(34) 
(35) 



where e is the electromagnetic coupling constant. The 
last relation H28|l follows from the facts that the con- 
densates (0 leave a U(1)q symmetry unbroken (that 

is; n^M77 = 0), and that the system is a Q-insulator 

in which Q-charged excitations are all gapped (that is; 

'^D,77 = . 

For later convenience we shall carefully look into how 
these analytical expressions (1^ and (PSll are structured 
from distinct contributions of particles and antiparti- 
cles. When quarks are massless, the quark propagator 
can be specifically separated into the particle and an- 
tiparticle parts (see Eq. (|40|) '). At the quark one- loop 



and 



['^A/]off(p-p) 



(36) 



D 

The 



from which we can easily make sure = 

as is well-known in the CFL phase l3^ 
particle-antiparticle contribution in the off-diagonal part 
[m^j]off(p-a) is vanishingly small, so we dropped it off from 
the above relations. 

It is of great importance to realize that the diagonal p- 
p contribution to the Meissner mass is negative one third 
of the Debye mass counterpart, while the p-a excitations 
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provide a positive contribution that is twice larger than 
the p-p's. Adding them together we finally acquire the 
Meissner mass contribution that is positive one third of 
the Debye mass contribution. In the off-diagonal part, in 
contrast, the situation is much simpler and the ratio is 
positive one third itself. 

In other words, in the diagonal part, the p-p loops al- 
ways tend to induce paramagnetism, while diamagnetism 
originates from the p-a loops. Usually in the super- 
conducting phase, the diamagnetic tendency is greater 
enough to bring about the Meissner effect. In gapless 
superconductors, however, antiparticles are never gap- 
less and only the p-p loops are abnormally enhanced due 
to gapless quarks and their large density of states near 
the Fermi surface causes the opposite phenomenon to the 
Meissner effect, i.e., chromomagnetic instability. Also, 
this can be seen from the relations (|34|l and (|36|l which 
hold in the presence of finite Ms or even in the gCFL 
phase. In the gCFL phase, a large density of states leads 
to a large Debye mass. Then, through (|34|l and H3()|) . it 
should be accompanied by a large negative contribution 
to the Meissner mass squared, which eventually results 
in chromomagnetic instability. 

As for (|^ . because the Ms dependence of antiparticle 
excitations is suppressed by Ms/fi, the relation is hardly 
changed for any Ms ^ fi. 



The quark propagator is divided into the particle and 
antiparticle parts by the energy projection operators 
Ap = i(l ± 7°7 • p). The exphcit propagator in the bd- 
gs sector, for instance, can be written down after some 
rearrangement of the Nambu-Gor'kov (1,2) (where 2 is 
assigned to the Nambu-Gor'kov doubler in our conven- 
tion) and color-flavor {bd,gs) indices, that is, in the basis 
(bd-l,gs-2,gs-l,bd-2), we have 



iS{p) 



A- C 
A+ C 
A 







A+ 
A, 












A+ 



p 


o" A: 



p 

A 







bd-gs 



gs-bd 



ip) 



sLgsip) 



^ gs-bd ip) 



(40) 



where the antiparticle part is 
1 



S\.b{p) = 



[po + ea-b) {po - £b-a) 

Po-{p + Hb) -iAi757° 
-i Ai757° Po + (p + ma) 



(41) 



II. COMPUTATION OF SELF-ENERGIES 

In this section we shall diagrammatically compute the 
calculable parts of one-loop self-energies for gluons and 
photon and derive the analytical expression for the sin- 
gular part of the Debye and Meissner masses. 

Using the mean-field quark propagator S{p), and the 
vertex matrices, F^, in color, flavor, spin, and Nambu- 
Gor'kov space, we can write the one-loop self-energies in 
a general form as 



1 



trT>iS{q+p)T^f,S(p). (37) 



The vertices F^ take a matrix form, 



-ig-f^T^S., 






igj^TaSi 



(38) 



for gluons {a — 1,...,8) where Ta are the SU(3)coior 
generators defined by the Gell-Mann matrices in color 
space and normalized as trT^T^ = ^Sap- For photon the 
vertex reads 
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-iejf'Q^ 






ieYQij 



(39) 



which is unity in color. We see that the flavor is not 
changed at any vertices but the color can be converted 
through the off-diagonal components in Ti, T2, T4, T5, 
Tg, and T7. We call the gluons corresponding to T3 and 
Tg as the color-diagonal gluons. 



and the particle part is 



Sa-b{p) 



1 



(po + £a-b) (po - £b-a) 

Po + (p-Ms) -iAi757° 
-«Ai757° Po~ (p- Ma) 



(42) 



In this representation the structure of the Dirac indices 
becomes much simpler and is given by A^ 7° and 757" 
attached with A. Then, it is important to notice that 
the (/X, ly) dependence in the integrand of (|37|l comes from 
the trace over the Dirac indices alone, that is 



diag(p- 
diag(p- 



.p)(p;q)=tr[7^A++p7VA+7"], (43) 
. .(p;q)=tr[7^A++p7VA;7°] (44) 



for p-p and p-a loops made of the Nambu-Gor'kov diago- 
nal components of the propagator H41() and (|42|l . For the 
Nambu-Gor'kov off-diagonal components we have like- 
wise; 



(p;g) =tr[7''A++p757''A+75] 



-'off(p-: 



(45) 
(46) 



A. Singularities in A\^; 



We will first discuss Ai and A2 that are degenerated. 
The generators in color, Ti and T2, have a non- vanishing 
component only between red and green and the flavor is 
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FIG. 6: An example of diagrams consisting of the diagonal 
components contributing to the self-energy for Ai,2- 



Ti T2 



ru 
/ 

\ 

gu 



gd 

N 

\ 

N 

/ 

rd 



FIG. 7: An example of diagrams consisting of the off-diagonal 
components contributing to the self-energy for Ai,2. Neither 
ru-gd nor gu-rd propagation involves gapless quarks, so that 
this diagram could not lead to any singular behavior in the 
gCFL phase. 



kept unchanged at the vertices, so that the self-energy 
for Ai^2 stems from four diagrams; 



diag(rg; s) 



diag(rg; u) 



diag(rg; d) 



off (ru-gd-.gu-rd 



(47) 



Here, diag(r5; s) means the quark loop composed of the 
diagonal component of rs and gs quarks (see Fig. ^ 
and dia,g{rg; u) and diag(rg; d) should be understood in 
the same manner. The last contribution, oS{ru-gd;gu- 
rd), represents the self-energy contribution from the off- 
diagonal components of ru-gd and gu-rd propagation, 
that is, the loop made up with ru quarks turning into 
gd and rd quarks turning into gu as shown in Fig. [T] 

In the case of Ai (and A2 equivalently) only ^'^^^^(^^g. ^-j 
may contain gapless quarks in the gCFL phase in both 
lines of the loop diagram, and is expected to provide a 
singular contribution. As a matter of fact, the singu- 
lar part presumably dominates over the screening mass 
behavior near the gapless onset as a function of Mg. 
This expectation is confirmed by the comparison between 
functional forms of only the rs-gs polarization and the 
full result derived from the potential curvature we are 
dealing with in the next section. From Fig. |S1 taking in 
the full results in advance, we can clearly observe that 
the physics near the gCFL onset is to be described by 
the rs-gs excitation only. 

The Debye screening mass, or n°°, has only the 
particle-particle part (and the antiparticle-antiparticle 
part that is negligible) because ^™g(p_a) (Pi 0) = 



'^'fi(p-a) (P' 0) = and is written as 



nd?ag(rg; s) il 



0) 



/ C^PP^'^d'iag(p-p)b;0)Z^diag(rg;s)(p) , 



(48) 
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FIG. 8: Comparison between the rs-gs contribution with an 
offset 0.39 and the full result for the Meissner mass squared 
for j4i,2. This means that the gross feature near the gCFL 
onset is provided by the rs-gs diagram alone. 



^hereTOOg(p_,)(p;0) 



2 and we defined 



^diag(rg:s) — /p-p [^bu-rs ! M, gS, 6m] 

+ fp.p[-~£rs-bu\bd,gs,bu] fp_p[ebd-gs;bu,rs,bd] (49) 
-I- /p_p [Sgs-bd; bu, rs, bd] . 

The p-p integrand /p_p is explicitly given by 

fp-p[eA-B;C, D, E] 



[ea-b+eb-a] [£a-b-£c-d] [£a-b+£d-c] 



tanh 



£a-b 



2T 
(50) 



with the energy dispersion relation £a-b defined in H15|l. 

In the CFL phase, as we have already mentioned, the 
bd-gs and rs-bu quark dispersion relations are identical 
and there are vanishing combinations in the denomina- 
tor of (inOJl like Egs-bdiq + p) - £rs-bu{p) ^ as q ^ 0. 
The numerator goes to zero in the same limit of g — > 
and consequently it amounts to the derivative with re- 
spect to e acting onto the distribution function, i.e., 
((9/9eA-_B) tanh[£^_B/2T]. At zero temperature this is 
proportional to the delta function d{eA-B) that takes a fi- 
nite value only when gapless quarks (that is; SA-Bip) = 
for some p) are present. 

Using the solution and H12|) and assuming Ai = 
A2 — A that is known to be a good approximation in 
the CFL phase, the polarization at zero temperature can 
simplify as 



ttOO 

-^Miag(rg; s) 



2e2 



A2 
2e3 



(51) 
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with p. = fi- M^/6n and e = \/ip~ P-)^ + A^. This 
expression is vahd only up to the onset where the gCFL 
starts. ActuaUy, right at the onset where Mf/2fi = A, 
the energy dispersion relation for gs-bd and rs-bu quarks 
is quadratic near the Fermi momentum and approxi- 
mated as e{p) - M2/2/i ~ (p - p-)^/2A. Then, the p- 
integration for the first term of (|51(l picks up the den- 
sity of states aX p = pL and ends up with an infrared 
divergence. We would emphasize that this singular be- 
havior originates from the electron pressure which en- 
sures /ig = 0. Otherwise, two quarks running through 
the loop diagram would not be gapless simultaneously at 
the gCFL onset. 

In the CFL region before the gCFL occurs, it is rather 
the second term of (|51ll which is a major contribution, 
which can be further approximated as 



ttOO 

diag(rg: s) 



5V 
47r2 



-—m„ 
2 3 



(52) 



with any logarithmic divergence neglected. The loga- 
rithmic terms are, in fact, given by ^ A^ ln[A/A] that 
is much smaller than the leading terms of order and 
corrections are ^1% for our choice of A = 800 MeV and 
Aq = 25 MeV. We will always drop any logarithmic 
terms ~ In [A/ A] throughout this paper. 

Hence, the contribution to the Debye mass in the CFL 
side is obtained as 



r 2 1 _ 2 

F7lDjdiag(r3; s) — ■ 



(53) 



Interestingly enough, this result exactly corresponds to 
one-flavor contribution (i.e. one-third) of l|30|) and does 
not have Ms-dependence until the gCFL onset. 

We can immediately extend our discussion to the 
Meissner mass once if we analyze the (/i,i^) structure in 
the polarization. The p-p contribution to the Meissner 
mass is 



*53)^diig(p-p)(p; 0) 



_ 1, 

3 ~ 3' 



(54) 



"'^'iag(p-p) (P! 0). 



In the third line we made use of averaging over the an- 
gle integration that leads to one third, that is, / d?'p{q- 
p)^U{p) = ^ Jd^pLi{p). In evaluating the Debye and 
Meissner masses Z^diag(rg; s) ip) is common, so that from 
the above relation we reach the general relation. 



rewrite /p_p into /p_a immediately, in view of the differ- 
ence between (|41|l and (I42|l . by replacing corresponding 
SA-B by EA-B, and p — /i^ by — (p+/i/i). A similar anal- 
ysis to (|54|) on the Dirac trace part leads to a factor as 
follows; 



lis. 



m3 



)^diig(p-a)(p;o) 



2[S'i ^fp>] ^l + {q-pf 



(56) 



Then, after some calculations along the same line as 
Refs. [23 I33 with the approximation that A and Ms 
are all neglected, we acquire 



j2 l(p-a) 



:— lim (5ii 

2 g^O^ 



diag(rg; s) 



(p-a) 



(57) 



1 2 A^ 
'2 "'9" 12^2 ' 



that is consistent with H35|l divided by the flavor factor 
A'f = 3 up to the ultraviolet divergence. 

After all, the sum over the p-p, p-a, and a-p parts 
finally amounts to 



[^2 1 _ r 2 i(p-P 

['''Mjdiag(rg; s) — I"' 



-I- 2 X ['nM]d'iag(rg;s) 



-^'^-ldiag(rg; s) 

(subtraction) 



(58) 



2 ^ 



in the CFL phase, where "subtraction" is a term added 
by hand to remove the A2-term which should be renor- 
malized. Eq. (|S7jl indicates 



A2 

(subtraction) — — - 
67r2 



(59) 



per one polarization diagram of Nambu-Gor'kov diago- 
nal particle-antiparticle excitations that should contain 
ultraviolet divergences. 



B. No singularities in ^4,5 and Ae,? 

The analyses in the previous subsection give us infor- 
mation about the origin of singularities around the gCFL 
onset. When two quarks involved in the loop diagram be- 
comes gapless at the same time, the divergent density of 
states compels us to face the imaginary Meissner mass. 
We can easily comprehend that this kind of singularity 
near the gapless onset is absent for ^4^5 and Aq,-;. 

The self-energies are diagrammatically decomposed as 



^2-1 (p-p) 
-/^'^ Jdiag(r(/; s) 



h 21 

2 V^D\d\a,g{rg; s)i 



(55) 



that is completely consistent with (|34|l . 

The p-a part generates ultraviolet divergent terms pro- 
portional to A2, A2ln[A/A], and (M2/^)2 ln[A/A] |33l |. 
In the present work, as in the treatment of the Debye 
mass, we ignore any logarithmic divergences. We can 



and 



1144 



aiag(or; 



n't." , 

diag(or; u) 



diag(fcr; d) 



ofi{hu-rs;ru-' 



diag(g6; s) 



diag(gb; u) 



diag(g6; d) 



ofi{gd-hs:hd-gs) ' 



(60) 



(61) 
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FIG. 9: An example of diagrams consisting of the diagonal 
components contributing to the self-energy for diagonal glu- 
ons and photon. 



bd 



bd 



FIG. 10: An example of diagrams consisting of the off- 
diagonal components contributing to the self-energy for di- 
agonal gluons and photon. 



None of these above polarizations is composed of would- 
be gapless quarks only, i.e. bd-gs and rs-bu quarks only. 
Therefore, we can conclude at this stage of analyses that 
^4^5 and Aqj do not exhibit chromomagnetic instability 
at least in the vicinity of the gCFL onset. 



C. Singularities in A^fi,- 



For color-diagonal gluons ^3, Ag; and photon A^, nei- 
ther color nor flavor is changed at the vertices. This 
means that 9 diagonal and 6 off-diagonal diagrams are 
possible. For example, the self-energy for ^g, for in- 
stance, has two types of singularities; one is associated 
with the bd-gs gapless quark dispersion relation. 



1 4 

bd-gs 2 diag(gg; s) g diag(fcb: d) 

2 

- - W" 

2 oii{bd-gs]bd-gs) ' 



(62) 



(64) 



L 88 J non-singular 3 diag(rr; 11) g diag(rr; d) 

114 

3 diag(gg; «) ^ 3 diag(£/g; d) ^ 3 diag(66; s) 
3 ofi{gu-rd\gu-rd) 3 o^{ru-gd-^ru-gd) 

2 2 

3 o^[gd-hs\gd-hs) 3 off (fcs-ru;6s-ru) ' 

The diagonal contributions, as shown in Fig. |51 are 
easily available from the previous results only by replac- 
ing the color-flavor indices appropriately and adding i to 
adjust the combinatorial factor. It should be noted that 
instability occurs in this case whenever any of quarks be- 
comes gapless because two virtually excited quarks are 
identical in the color-diagonal channels. 

The off-diagonal expression is also deduced simply by 
means of proper arrangement of the color-flavor indices 
and replacement of the numerator of H5U|I as ea-b i (p ~ 
^lb) ^ —i^AB- Because ea-b ± {p ~ I^b) — ^ab in the 
vicinity of the Fermi surface that is dominating over 
the p-integration, the functional form of the off-diagonal 
contribution is essentially the same as the diagonal con- 
tributions. There are two important features to remark 
here, though. 

First, the self-energy contributions coming from the 
off-diagonal loop are free from ultraviolet quadratic di- 
vergence even in p-a excitations. It is because the off- 
diagonal propagator is proportional to, not ea-b ± (p — 
/is), but IS.AB that is not rising with increasing p. Thus, 
the subtraction term proportional to A^ is no longer nec- 
essary. 

Second, from the trace over the Dirac indices H45|l and 
(|45|l we notice that 



X. 



off 



-'off ~ -'diag 



(65) 
(66) 



for both p-p and p-a cases, that means the diagonal and 
off-diagonal diagrams contribute to the Debye and Meiss- 
ner masses differently. Using (|65ll , (|66f) and from a similar 
analysis to (|54|l . we immediately see the relation. 



and the other is associated with the rs-bu gapless quark 
(almost quadratic) dispersion relation. 



1 4 

rs-bu 3 diag(rr; 5) 3 diag(b6; u) 

2 

2 ofi{rs-bu;rs-hu) 



(63) 



The one-loop diagrams corresponding to H^j'^g^fif,. ^-j and 



oif{bd-gs;bd-gs) 



are shown in Figs.l^andllUlto take some 

examples. 

For completeness we shall enumerate the rest here. 
The non-singular contributions to the Ag polarization are 



i^l-AoSibd-gs-M-gs) — g ['7l|)]off( bd-gs ;6d-gs) > (67) 

that is consistent with H36|) . The sign difference between 
(|55|l and 1)6 7|) originates from the negative sign in H65|l . 
The p-a off-diagonal contribution to the Meissner mass 
squared turns out to be negligible. 

In addition, the negative sign difference between H65|) 
and H66|) can explain the results for the "neutral" and 
"charged" condensates discussed in a two-species model 
in Ref. j^^. As we have mentioned above, the func- 
tional form of IIj." , w n'f.'^ ,w and II'^" 

J f — ^.y diag(bb; d) ' 



diag(gg; s) ' 



ofi{bd-gs:bd-gs) 



are almost identical up to the combinatorial factor i. 
Let's consider the "neutral" case fictitiously following 
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the argument in Ref. |25j, for which we assume that 
bd and gs quarks had the opposite charge correspond- 
ing to the two-species model composed of bd and gs. 
Then, the polarization related to the screening masses 
would be n'*" — H^'^ + 11'^'^ — 11'^'' 

diag(5t/; s) diag(fcb; d) off(bd-gs-,hd-gs) ' 

From (|65|) we see that the Debye mass positively di- 
verges at the gapless onset, while the singularities in 
the Meissner mass cancel to vanish. In the "charged" 
case, on the other hand, we assume that bd and gs 
quarks had the same charge. In the same way, the 
polarization leading to the screening masses would be 

n^'' = m.'' , , + n'f.'' , -I- n^^,, , , , , . obvi- 

diag{gg: sj aiag(oo; a) on{ba-gs:Da-gsj 

ously, the Debye mass has no singularity at the gapless 
onset, but the Meissner mass negatively diverges. This 
is exactly what was found in Ref. In QCD, how- 

ever, the diquark condensate is actually a mixture of the 
"neutral" and "charged" condensates and the Debye and 
Meissner masses both diverge positively and negatively 
from the common origin at the gCFL onset. 



As a final remark of this section, we shall draw at- 
tention to the contrasting difference between the singu- 
lar parts (|62|l related to bd-gs quarks and H()3I) related 
to rs-bu quarks. In general, the Meissner mass nega- 
tively diverges on the bd-gs, rs-bu, and gu-rd boundary 
lines drawn in Fig. 01 When Ms grows larger, as seen 
in Fig. |S1 the blocking region of the bd-gs dispersion re- 
lation becomes wider and the singular character in H62() 
weakens. In contrast, fj,,, is kept to be very close to the 
rs-bu line in Fig. 21 for all Mg and in all the gCFL region 
(|63|l remains to provide a negatively huge contribution to 
the Meissner mass. If it were not for electrons, fie would 
be within the bandgap and ]<j'3\i would not be singular 
at all. In this sense, the singularity of Ht)3|l is induced 
by the presence of the electron terms which force /ig to 
stay in the singular region slightly below the dashed line 
in Fig. 01 To evaluate that huge contribution correctly 
from it is essential to resolve a tiny blocking re- 

gion Sp ^ /ig/3/i^ with a great accuracy, which is next to 
impossible technically. 



Although the strength of singularities 1)62(1 and Ht)3|) 
at the gCFL onset is exactly equal at T = 0, the tem- 
perature dependence is drastically different; as seen from 
Fig. ^2 Ht)3|) is no longer singular even at T = 0.5 MeV. 
This is because at finite temperature rs-bu quarks not 
from the blocking region but from thermal excitations 
can cancel the Q-charge driven by electrons. There 
needs not to be a blocking region associated with rs- 
bu quarks therefore. Actually, the rs-bu quark density is 
estimated as V^nATp.'^ /tt"^ when its dispersion relation 
is quadratic, and it can balance out the electron density 
/ig/37r^ even at T = (/ig/3/i)^/27rA that is as tiny as ~ 
eV for /ie ~ a few MeV. Consequently there is no gapless 
rs-bu quarks any more already at tiny temperatures. 
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FIG. 11: Each contribution from possible singular diagrams 
at T = and 0.5 MeV. In the figure hd-gs and rs-bu repre- 
sent contributions to the Meissner mass from [ngg]t,d_gs and 

L 88 jfs-bu- 



III. POTENTIAL CURVATURE 

In order to derive the information on the full evalu- 
ation of screening masses, let us elaborate a numerical 
method based on the previous analytical consideration 
to compute the Debye and Meissner masses here. It is 
well-known that the Debye mass can be expressed as the 
potential curvature as follows [s^l ; 



I'D, 0/3 



(68) 



where /x's are the color chemical potential coupled to the 
generators Tq,. That means, /it's here involve off-diagonal 
components in color-flavor space besides The ther- 
modynamic potential defined with all //'s is denoted by 
here. The potential curvature ((68|l is evaluated at 
vanishing /i's except and /ig which are fixed by the 
color neutrality conditions. In QCD at finite tempera- 
ture, moreover, pure gluonic loops produce a Debye mass 
^ gT , which is not included in Eq. H68|l. 

All we have to know to evaluate the thermodynamic 
potential are the energy dispersion relations that are 
much easier than calculating the loop integral (|37() di- 
rectly. We show the Debye masses as a function of //z 
inferred from ((68|l in Figs. 1121 and 1131 The singular be- 
havior associated with the gapless onset emerges in an 
upward direction, that is accounted for from the rela- 
tion 134(1 . In other words, through the relation ((34(1 . the 
positively rising results for the Debye masses squared im- 
ply negative Meissner masses squared for Ai^2 and ^3,8,7 
around the gCFL onset and for all the gluons for large 
M'^/fj,. It is worth noting that the Q-photon has a finite 
Debye mass for large M^/fJ. because the system in the 
gCFL phase has a nonzero density of electrons that can 
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FIG. 12: Debye mass for Ai^2 (solid), A^^^ (dashed), and Ag,? 
(dot-dashed) obtained from the potential curvature at T = 0. 



terms diverging like ~ , and VIa is the thermodynamic 
potential defined in the presence of the gauge fields Ai 
(Note that bold symbols stand for matrices in color-flavor 
space in our notation.) in the Lagrangian ^ as 



(70) 



Equation H69(l needs some more explanation. It should 
be worth mentioning that H69I) is a well-known general 
relation if we do accomplish the angle integration with 
respect to the loop momentum p. However, this is tech- 
nically difficult, for p makes an angle with on evalu- 
ating Ha- The benefit of Ht)9|l is that we can reduce the 
p-integration into a one-dimensional pa-integration as if 
the rotational symmetry were not affected by Ai. We 
shall take advantage of this simplification at the price of 
treating Ms as chemical potential shifts. 

This relation can be proved in the following way. When 
the momentum p in the loop-integral is chosen to be in 
the 3-direction, we can show that. 
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FIG. 13: Debye mass for the color-diagonal gluons A3, Ag 
and photon A-y at T = 0.5 MeV. Three solid curves represent 
three eigenvalues, and dotted curves are the Debye masses 
with mixing turned off by hand. The electromagnetic coupling 
constant is chosen to hee/g = 1/2, which is presumably larger 
than a realistic value, to make visible the mixing effect of A^. 



screen Q-charge. 

Once we develop a relation connecting the thermody- 
namic potential to the Meissner mass like the above rela- 
tion (|B8|I with respect to the Debye mass, we can numer- 
ically compute the Meissner mass including not only the 
singular parts but also all the contributions. As a matter 
of fact, we have found a relation. 



('™M,a/3)ba 



1 ^ 



(69) 



A=0 



where "bare" means the Meissner mass squared with 
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(71) 



indicating that H'^'^ contains only the p-p contribution 
and and 11^^ the p-a contribution. From these ex- 
pressions in the same way as (|54f) we eventually arrive at 
the relation. 



Mn"/3(0'0) 



n^^^(o,o) 



n33,(0,0)]. (72) 



We should remark that, even though we do not know 
an explicit expression for the ru-gd-bs quark propagator, 
the {fi, v) structure comes out only through the energy 
projection operators and the above argument is gener- 
ally valid for any quark propagator as long as can 
be treated as an effective chemical potential shift. The 
relation (|72|l is readily rewritten in terms of the thermo- 
dynamic potential, which leads us to H()9|l . 

In our model study, the four-fermion interaction is non- 
renormalizable and ^a suffers ultraviolet divergent terms 
proportional to A^ . Thus the Meissner masses derived di- 
rectly from the potential curvature should have those di- 
vergences to be subtracted properly by hand. As argued 
in the previous analytical study, we conjecture that the 
subtraction term is simply inferred by counting divergent 
diagrams, which results in; 



''Ma 



(m 



M.QQ jbarc 



3 X 



— ('^M,77)bare + 4 X 



4/,a/3 = (»7tif,a/3)barc , (1 < a 7^ /3 < 



"6^ 



67r2 



(1 < a < 8) (73) 
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FIG. 14: Meissner mass squared for Ai^2 as a function of 
Ma/^ at the temperatures T = (solid), 0.5 MeV (dashed), 
and 5 MeV (dot-dashed). The vertical dotted line represents 
the gCFL onset. The model parameters are chosen to yield 
Ao = 25 MeV at zero temperature and zero strange quark 
mass. The solution of the coupled equations ceases to exist 
at Wi/^ = 144.7 MeV and 140.5 MeV for T = and 0.5 MeV 
respectively. 



These are all determined simply from the diagrammatic 
counting of the Nambu-Gor'kov diagonal p-a loops and 
the associated subtraction H59|) with proper coefficients. 
In the case of for example, the divergence factor in 
unit of A^/67r^ coming from is ^xlj + f) — | where 
^ is a combinatorial factor. In the same way and 
()64|l give I and | respectively. The sum of these results 
in 3 X (A^/67r^) as in l(7^ . It is interesting to see that the 
divergent terms cancel each other in the color-mixing and 
photon-mixing channels. From these above expressions 
we have obtained the results as presented in Figs. 11411151 
CSl and El 

Now, before addressing the instability problem, let us 
compare our numerical results shown in Figs. 1141 1151 1161 
and 1171 to the known analytical estimates. In view of 
our results at = 0, all the gluons except Ag have 
the identical Meissner mass 0.892TOg which is consistent 



with (|25|l . We have confirmed that the agreement be- 
comes better for larger /j, ^ A because the analytical es- 
timates require this. The mixing channel with e/g — 1/2 
(chosen presumably larger than realistic values to make 
the mixing effects apparently visible) results in 1.1 Sto^ 
for As and zero for A^, which is also consistent with the 
analytic estimate m^ss ~ 1.1 5m^ from H27|) . All these 
suggest that our method works well. 

Figure H"!] shows the Meissner mass squared for Ai^2 as 
a function of / fj,. The vertical dotted line corresponds 
to the gCFL onset and is located near the onset condition 
given kinematically (M^ / /Lt)c ~ 2 A. (We chose the model 
parameters to yield Aq = 25 MeV for = T = 0.) We 



FIG. 15: Meissner mass squared for A4^5 as a function of 
M^/fi at r = (solid), 5 MeV (dashed), and 10 MeV (dot- 
dashed) . 
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FIG. 16: Meissner mass squared for Ae,? as a function of 
M^/fi at r = (solid), 5 MeV (dashed),' and 10 MeV (dot- 
dashed) . 



see that the Meissner mass squared is negative in the en- 
tire gCFL region at zero temperature, that means chro- 
momagnetic instability. The zero temperature behavior 
for Ai,2 is consistent with the results reported in Ref. 
except for the asymptotic property that the Meissner 
mass in Fig. El does not approach zero for large /fj, 
as it does in Ref. . The singularity around the gCFL 
onset is drastically smeared by the temperature effect as 
seen in the curve for T = 0.5 MeV, though chromomag- 
netic instability stiU exists for Mf/fi > 84 MeV. At 
higher temperatures, as shown by the T — b MeV curve, 
the Meissner screening mass smoothly goes to zero with 
increasing M'^/fi since the system then comes to reach 
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FIG. 17: Meissner mass squared for three eigenmodes com- 
posed of A'i, An, and A-i as a function of Ml j ^ at T — 
0.5 MeV. The electromagnetic coupling constant is chosen 
to be e/g — 1/2. The dotted curves are for ^3, As, and A.y 
with mixing among them turned off. 



the normal phase without any first-order phase transi- 
tion. (The T = 5 MeV results in fact have a tiny unstable 
region at M^/n = 125 MeV, see Fig. [2) 

In Fig. (and in Fig. I16|l we plot the results for ^4^5 
(and for Aqj respectively). The gross features are in 
good agreement with the results presented in Ref. ^ 
and all these gluons do not suffer chromomagnetic insta- 
bility until M^/fj, > 105 MeV at zero temperature. We 
remark that the instabilities for //i > 100 MeV persist 
even at T = 5 MeV, but it disappears at T = 10 MeV. 

The Meissner masses corresponding to A3, As, and 
A^, as we have exposited before, are negatively diver- 
gent if any of quark excitation energies take the quadratic 
form as a function of the momentum. This situation 
takes place only at the gCFL onset. The results inside 
the gCFL regime are negatively large but not divergent 
actually. In fact a finite density of electrons requires 
an almost but not exactly quadratic dispersion relation 
anywhere in the gCFL phase. The outputs accordingly 
have severe dependence on the precise form of the almost 
quadratic dispersion relation, that is, the deviation from 
the quadratic form, which is hard to evaluate numerically 
with an extreme accuracy. We have confirmed the above 
mentioned behavior quantitatively in our calculations, 
that is, negatively large Meissner masses squared for ^3, 
As, and A^ in the entire gCFL region, but we could not 
reduce tremendous numerical ambiguity. Instead of tack- 
ling numerical difhculties that would be milder at finite 
temperature, we have calculated the Meissner masses at 
T = 0.5 MeV. In Fig.Elthe dotted curves are the masses 
without mixing taken into account and the solid curves 
are the eigenvalues of the 3x3 mass squared matrix. 
We clearly see that instability certainly occurs in one or 
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FIG. 18: The correspondence between chromomagnetic insta- 
bility for ^4,5) and Aej and the CFL, uSC, and 2SC phases 
as well as unpaired quark matter (UQM). The temperature is 
chosen to be T = 5 MeV to draw this figure. 



two of three eigenmodes. Furthermore, we note that any 
effect of the almost quadratic dispersion relation cannot 
be perceived from the T = 0.5 MeV results because it is 
already smooth enough at that temperature as is obvious 
also in Fig. ^] 

The unstable regions read from these figures are to be 
considered as slices with T fixed on the phase diagram. It 
is instructive to compile all these slices and map respec- 
tive unstable regions onto the phase diagram. Figure [T51 
is one example to locate the phase boundary and the un- 
stable region and to see them on a common basis. In 
this way we clarified the unstable regions on the phase 
diagram as we have already discussed at the beginning 
of this paper. 



IV. SUMMARY AND OPEN QUESTIONS 

In summary, in dense neutral three-flavor quark mat- 
ter, we figured out the computational procedure to de- 
rive the Debye and Meissner screening masses from cur- 
vatures of the thermodynamic potential with respect to 
gluon source fields. We investigated chromomagnetic in- 
stability by changing the temperature T and the Fermi 
momentum mismatch M'^/fi and explored the unstable 
regions for respective gluons and photon on the phase 
diagram involving the (g)CFL, uSC, and (g)2SC phases. 

We can say that the instability lines we added on the 
renewed phase diagram tell us, at least, at what /i or 
Ms the homogeneous (g)CFL, uSC, and (g)2SC phases 
remain to be stable at a certain T. From the present 
work, however, we cannot extract any clue about what is 
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FIG. 19: The unstable region in which at least one gluon has 
imaginary Meissner mass is shown enclosed with solid thick 
lines, which can be interpreted as boundaries between the 
homogeneous and (colored) crystalline phases. 

going on inside the unstable regions. One possible solu- 
tion would be a crystalline-like phase as is suggested in 
the two-flavor case . If this is true also in three-flavor 
quark matter, it is most likely that the instability lines 
are to be regarded as the boundaries that separate the 
homogeneous and crystalline superconducting phases. In 
this sense, the boundaries drawn by thick lines in Figs.^ 
121 and|3have a definite physical meaning even after a true 
ground state inside the unstable regions will be identi- 
fied. Strictly speaking, once chromomagnetic instability 
occurs in any gluon channel, it would affect the other 
gluon channels. Consequently instability lines with re- 
spect to one gluon going inside the unstable region with 
respect to another gluon might be substantially different 
from what we have seen in this work. They are actually 
the most outer boundaries only that separate the homo- 
geneous and crystalline phases, as shown in Fig. 1191 

Our approach here is essentially along the same line 
as in Ref. and so this work leaves the same open 
questions enumerated in Ref. how gauge field fluc- 
tuations affect the critical point and the order of re- 
spective phase transitions what the nature of K'^- 
condensation in the gCFL phase is and how it affects 
instability 36, 37], what difference the 't Hooft (instan- 
ton) interaction makes that induces an interaction like 



jAapMs, where the chiral phase transition is located 
and how our (M'^ / fj,-T) phase diagram is mapped onto 
the (/i-T) plane with Ms solved dynamically. 

In addition to these above mentioned issues, we should 
make an improvement how to take account of Mg fully. 
The derivation of our formula needs the energy pro- 
jection operator for massless quarks and works out as 
long as the Mg effects can be well- approximated by a 
chemical potential shift. From this reason, we could 
not augment the phase diagrams for stronger couplings 
(Figs. 16 and 17 in Ref. ^ for example) with the insta- 
bility lines. The improvement could be implemented by 
the full angle-integration in p to acquire ft a ■ 

From the theoretical point of view, logarithmic diver- 
gences that have been simply neglected in the present 
work must be taken seriously. In particular, logarithmic 
terms depending on Ms seem to remain even in the nor- 
mal phase ll^, while in our work the Meissner masses 
turn out to be zero properly in the normal phase, sug- 
gesting some cancellation mechanism between Ms depen- 
dent divergences. This must be further investigated in 
a field-theoretical way. Another problem of theoretical 
interest is how much singularities are inhibited by non- 
perturbative resummation. However small the coupling 
constant g is, the polarization is divergently large at the 
one-loop level, and generally, that demands a sort of re- 
summation, for instance, by means of the ladder approx- 
imation. The instability might not be cured itself by 
resummation, but the singularities near the gCFL onset 
would be possibly smoothened. 

Phenomenologically, since not only the gCFL but also 
the g2SC phase is within our scope especially once we 
will be able to approach the stronger coupling case, it is 
inevitable to examine the competition between the mixed 
and crystalline phases around the g2SC region. The en- 
tanglement between these phase possibilities should play 
an important role to complement the phase diagram for 
the large Ms or moderate density region [38l |. 
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